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Abstract. The angular correlations measured in proton-proton collisions at y^s = 7 TeV are decomposed 
into contributions from back to back emission and elliptic flow. Modeling the dominant term in the corre- 
lation functions as a momentum conservation effect or as an effect of the initial transverse velocity of the 
source, the remaining elliptic flow component can be estimated. The elliptic flow coefficient extracted from 
the CMS Collaboration data is 0.04 — 0.08. No additional small-angle, ridge-like correlations are needed to 
explain the experimental data. 



' 1 Introduction 

, Proton-proton collisions at the CERN Large Hadron Col- 

■ lider with the energy of 7 TeV are the most violent ele- 
' mentary collisions studied in a laboratory. The high mul- 
. tiplicity of produced particles means that a short-lived, 
' high density system is formed. By analogy to the physics 
, of ultrarelativistic heavy-ion collisions, one could expect 

■ new collective phenomena to emerge in the highest mul- 
] tiplicity events. A striking experimental observation has 
, been made by the CMS Collaboration, measuring two- 
' particle correlations R{A(j), Arj) in the relative azimuthal 
, angle and the relative pseudorapidity [1]. In high multi- 

■ plicity events an enhancement of the correlation on the 
I near side (Z\0 ~ 0) is observed, extending in pseudorapid- 
. ity in the form of a ridge. This gives rise to speculations 
' on the origin of the effect [2,3]. 

We present a quantitative analysis of the angular cor- 
relation R{A(j)) measured by the CMS experiment for pair 
of particles separated in pseudorapidity |Z\77| > 2.0. The 
angular correlation function is decomposed into a contri- 
bution from global momentum flow and a smaller elliptic 
component. The dominant structure in the correlations 
function in the relative azimuthal angle is due to an en- 
hancement of back to back emission. This effect is qualita- 
tively described by Monte Carlo event generators [1]. The 
experimental data compared to the PYTHIA generator re- 
sults show two essential differences. For all transverse mo- 
mentum and multiplicity classes presented, the measured 
angular correlations show deviations from model calcula- 
tions. In simulations and in the data, the correlations are 
enhanced for the emission with A(j) — tt, but the width 
and the height of the peak is different. The most impor- 
tant new feature in the measured correlation functions is 
the appearance of a small-angle enhancement of the dis- 
tribution. This small-angle structure (ridge) is present in 
a broad range in the relative pseudorapidity between the 



particles in the pair. This small-angle enhancement is not 
due to the usual correlations from the jet fragmentation, 
and is not reproduced by any Monte Carlo event genera- 
tor. 

We analyze the possibility that the new effect is due 
to the elliptic flow present besides the dominant back to 
back angular correlations. For some transverse momen- 
tum and multiplicity cuts, the relative importance of the 
elliptic flow component is big enough to make it explic- 
itly visible in the small-angle region. However, the second 
harmonic cos(2Z\0) cannot be unambiguously unraveled 
in the whole range (O-tt) from the dominant background 
in the correlation function. 

The importance of the elliptic flow is that it is a signa- 
ture of the collective expansion of the fireball created in 
relativistic heavy- ion collisions [4] . It is important to check 
whether similar phenomena occur in elementary collisions, 
but the search for the elliptic flow in smaller systems is 
hindered by statistical fluctuations and non-flow correla- 
tions. Several model estimates of the elliptic flow gener- 
ated in proton-proton collisions have been given [5-13]. 
Most commonly, one assumes that the elliptic flow is gen- 
erated during a short hydrodynamic expansion, in a very 
similar way as in heavy-ion collisions. The most impor- 
tant differences between calculations concern the origin 
of the azimuthal asymmetry of the initial energy density 
distribution in the transverse plane. If the presence of an 
elliptic flow of collective origin could be demonstrated in 
proton-proton collisions, it would indicate the a strongly 
interacting flreball has been formed in hadron collisions at 
the highest energies. 

The elliptic flow correlations are subleading. The ex- 
traction of the elliptic flow coefficient requires a careful 
analysis of other effects. In the following, we propose two 
models for the dominant back to back correlations. The 
first estimate is based on the momentum conservation ef- 
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fects, requiring the conservation of the total transverse 
momentum. The second calculation assumes that particles 
separated in pseudorapidity are emitted from sources mov- 
ing (on average) in opposite transverse directions. This 
picture is natural if the two particles come from two dif- 
ferent back to back jets. In a high multiplicity event, and 
for low transverse momenta of the particles, a simple pic- 
ture can be used, with two sources moving with opposite 
transverse velocities. 

To reproduce the observed ridge correlations, an addi- 
tional azimuthal asymmetric component (elliptic flow) is 
added. The estimated elliptic flow coefficient is t'2 = 0.06- 
0.1, for the parameters where the ridge is observed. The el- 
liptic flow interpretation agrees qualitatively with the mul- 
tiplicity and transverse momentum dependence of the ob- 
served effect. The strength of the small-angle two-particle 
correlations is proportional to V2, and should increase as 

p1- 



2 Momentum conservation 

The correlation function R{A(j)) is obtained selecting par- 
ticles separated in pseudorapidity [1] . This procedure elim- 
inates jet-like correlations on the near side. The form of 
the dominant remaining correlations shows, that the as- 
sociated particle is emitted preferentially in the opposite 
direction in the azimuthal angle. This could be a conse- 
quence of the momentum conservation in the microscopic 
particle production mechanism. Such correlations can be 
very important in small multiplicity events, as in proton- 
proton collisions [14-17]. In the following we take into ac- 
count the transverse momentum conservation in the same 
way as proposed in [14, 17]. In an event of multiplicity M, 
the multiparticlc! distribution in transverse momenta Pi 
and pseudorapidities r}i is given as 



/m(Pi,??i, • • --.Pm^Vm) = 
S'^{pi + . . .Pm)/(pi, ??i) • • • f{PM, Vm)/ 



6'^{pi + . ..pm) 



/(pi,77i) . . . fipM,riM)d^Pidrii . ..(fpMdrjM (1) 



The distribution / is 



.nP,V)^fav{p±,V){^ + 2v2 COS(20)) 



(2) 



and is normalized to one. M is not necessarily the total 
multiplicity in the event, it is rather the number of parti- 
cles created in the microscopic process conserving the mo- 
mentum. V2 is the elliptic flow coefficient, that is assumed 
to have negligible variation over the transverse momen- 
tum, pseudorapidity or multiplicity bins where the aver- 
ages are taken. The two-particle distribution is obtained 
from Eq. (1) by integrating over the momenta of M — 2 
particles. For large M, small V2, and to the order 1/A^, 
one obtains [14, 17] 



h{Pi,m,P2,ri2) = f{.Pi,Vi)f(.P2,V2) 1 + 77 



pI 2piP2C0s((/)i - (/>2)' 



M{pI)f 



(3) 



{■ ■ ■)f denotes the average over the full phase space. The 
single particle distribution requires one more integration 
over p2 and 772 



A(P,.) = /(p,.)(i + ^-^) . 



(4) 



The correlation function R{A(j)) measured by the CMS 
Collaboration can be expressed as 

R{A(P) = {<N> -1) 

/g d'^pidl]i(fp2d'n2f2{Pl,Vl,P2,V2)S{A(j) - 01 + 02) 

(5) 



Jg d'^pi drji d'^p2dr]2 fi {pi,T]i)fi{P2,V2)H^'P- 

-1) , 



where the integration is performed over the considered 
transverse momentum bin and fulfilling the condition on 
the pseudorapidity separation 2.0 < \r]i — 772] < 4.8, also 
the numerator and the denominator in the above expres- 
sion are averaged over events in a given multiplicity (A^trfe) 
class. We obtain a simple expression 



R{A(I)) = -ci cos{A(j)) + C2 cos{2A(j)) 



where the coefficients are 



Cl 



;i 111 1 



{pDpNeff 



C2 = 2((iV) 



eff 



l)vl 



(6) 



(7) 



(8) 



To extract the elliptic flow coefRcient the average multi- 
plicities {N) in each bin are used. These multiplici- 
ties are obtained from the integration of the p±_ distribu- 
tion of particles produced in proton-proton collisions at 
= 7 TeV [18]. For minimum bias events the transverse 
momentum distribution is well described using a Tsallis 
distribution 



dN 



dr]p±dp_i 



1 



nT 



(9) 



M M{j)\)f 



where the slope parameter T = 0.145 GcV and n = 6.6. 
For our estimate, we assume the same distribution in the 
whole pseudorapidity range, with the constant C adjusted 
to reproduce the total multiplicity in the interval jryj < 2.4. 
A noticeable increase of the mean transverse momentum 
with multiplicity is observed [19]. To take this effect into 
account, we assume that the slope parameter depends 
on the multiplicity Ntrk in the CMS acceptance region 
(|r;| < 2.4, p± > 0.4 GeV) T = (0.02 + My/N^) GeV. 
Such a distribution, when recalculated for the ATLAS ac- 
ceptance, reproduces the increase of the mean transverse 
momentum with multiplicity [19]. Using the parameteriza- 
tion (9) we calculate the mean multiplicities in different p± 
bins and multiplicity classes Ntrk- The events used in the 
analysis have at least two particles in the chosen kinematic 
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Fig. 1. The correlation function R{A(I)) as function of the relative angle (5), for particle pairs with pseudorapidity separation 
2.0 < |Zi??| < 4.8, for different multiplicity and transverse momentum bins [1]. The solid lines represent the fit — ci cos(zl(/)) + 
C2 cos(2Zi(/>), the dashed lines represent the — Ci cos(Zi(/>) term only. 



range. For small multiplicities we take the mean multiplic- 
ity (N) only for events fulfilling the condition that N > 2. 
The correction is made using a Poisson distribution. 

The expression (6) is used to fit the angular correla- 
tion function R{A(j)) measured in proton-proton collisions 
at = 7 TeV [1]. The quality of the fit is excellent. The 
angular correlation has two components, the momentum 
conservation and the elliptic aziniuthal asymmetry corre- 
lations. These effects account for all the angular structure 
observed in the experiment for particles emitted in differ- 
ent pseudorapidity regions. Some deviations are visible in 
the highest momentum bin, especially in the low multi- 
plicity class. This may signal the presence of two-particle 
correlation not of the assumed form (6). The expansion 
used in deriving momentum conservation effects breaks 
down for large p± and small Ntrk {M). The coefficient 
ci (~ 0.4-0.9), denoting the strength of the momentum 
conservation correlations, is of the order one, as expected 
from the formula (7). 

The extracted C2 can be used to calculate the strength 
of the elliptic flow in the particle azimuthal distributions 



V2 



C2 



2((iV) - 1) 



(10) 



which is equivalent to the second order cumulant method 
[20]. In Fig 2 is shown the elliptic flow as function of the 
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Fig. 2. The elliptic fiow coefficient «2 as function of the trans- 
verse momentum (Eq. 10), for the four multiplicity classes con- 
sidered. The lower open squares represent the elliptic flow coef- 
ficient extracted using the small angle data Acj) < ^ , assuming 
a constant background (Eq. 11). 



transverse momentum bin for the four multiplicity classes 
considered. The first observation is that the results for low 
multiplicity events differ from the others. It may indicate 
that the correlations in low multiplicity events are of a 
different origin than in high multiplicity classes. For the 
high multiplicity events, the value of the elliptic flow in- 
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creases from 0.03, in the lowest, to 0.1, in the highest p± 
bin. This value is too high to be interpreted as entirely 
due to collective effects. The strong increase with p±, and 
the large absolute value indicate that at least part of the 
effect is due to the non-flow correlations that cannot be 
separated. Hydrodynamic model estimates for the elliptic 
flow coefficient in proton-proton collisions are of the order 
0.04-0.06, if viscosity effects reducing the anisotropy are 
taken into account [5,13]. 

We have followed the same procedure, using a fit with 
the first and second harmonics, for the PYTHIA simu- 
lated distributions [1]. The resulting elliptic flow coeffi- 
cient is of the same order as found in the data. This indi- 
cates that the non-flow, jet-like correlations, included in 
PYTHIA have both cos(Z\0) and cos(2Z\0) components. 
Of course, the elliptic flow component in PYTHIA is not 
of a collective origin. On the other hand, the ridge struc- 
ture is not present in the correlation function from the 
Monte Carlo generated events. Moreover, since the parti- 
cle pair is separated in pseudorapidity, most of the short 
range correlations, coming from jets or resonance decays, 
are strongly reduced [21]. The remaining non-flow cor- 
relations are modifying the angular correlation function 
for large relative angles, as seen in the PYTHIA simu- 
lations. It indicates that the assumed form of the back- 
ground — ci cos{A(j)) is oversimplified, and that non-flow 
higher harmonics are present in the observed data. The 
small-angle structure can be explained as resulting from 
the elliptic flow, but the extracted values of the V2 coef- 
ficient depend on the form of the correlation R{A<j)) in 
the whole range O-tt, and contain a large contribution of 
non-flow origin. To estimate the value of the elliptic flow 
related to the appearance of the ridge-like structure, we 
fit a simple formula 



i?(Z\(/)) = b + C2 cos(Z\0) 



(11) 



in the range Z\(/) < The results for the multiplicity 
class 90 < Ntrk < 110 is shown as the lower points in Fig. 
2. 

The momentum resolution of the correlation data is 
not sufficient to look for a hydrodynamic origin of the 
effect. Future data with more statistics could allow to ex- 
tract a detailed behavior in the soft momentum region 
p± < 2.0 GeV. The charge independence of the effect 
found in [1] is consistent with a collective flow origin of 
the correlations. The multiplicity dependence of the sec- 
ond harmonic coefficient in the correlation function R(A(f)) 
(c2 c>c (N)) indicates that the underlying correlations {V2) 
are of a non-flow origin. On the contrary, the approximate 
multiplicity independence of the coefficient ci shows that 
the underlying correlations decrease linearly with the size 
of the system, which suggests a non-flow origin of the ef- 
fect. 

Let us comment on the first harmonic term in the cor- 
relations function R(A<j>). The formula (8) allows to esti- 
mate the value of effective multiplicity N^f / of the process 
for which the momentum conservation occurs. We find val- 
ues between 10, in lowp± bins, and 80, in high momentum 
and high multiplicity bins. Such values are small compared 



to the multiplicity in the range \ri\ < 2.4, and even more so 
with respect to the total multiplicity in the event. This, 
and the fact that higher harmonics in the angular cor- 
relation function are important shows that the observed 
dominant correlations have a more complicated form than 
cos(Z\^), especially in high the p± bins. 



3 Transverse source velocity 

In this section we analyze the angular correlations assum- 
ing a different model for the dominant back to back peak. 
In the first collision between partons a hard momentum 
transfer occurs. Subsequent particle production is deter- 
mined by the initial transverse momentum exchange. We 
use a simple estimate for the angular correlation resulting 
from such a mechanism. Particles separated in pseudora- 
pidity originate predominantly from initial partons with 
opposite transverse momenta. Let us assume that in the 
rest frame of the source connected to a given fragmenting 
parton, massless particles are emitted isotropically. If the 
source is moving with velocity f3± in the direction (p = tt 
the final angular distribution takes the form 



dN 



(12) 



d4> 1 + /3_L cos((/)) 
The source from which the other particle in the pair orig 




Fig. 3. Schematic picture of the source created in a proton- 
proton collision. At each rapidity it has a transverse velocity. 
The source is azimuthally deformed and its collective expansion 
would lead to a nonzero elliptic flow. 



inates, moves on the average with velocity — (Fig. 3). 
If the transverse momentum and pseudorapidity depen- 
dence factorizes from the angular one, we obtain for the 
angular correlation of the pair 



R{A^) = {{N) - 1) 



2 - /?2 + /3_L C0S(Z\(/)) 



(13) 



The above correlation function is peaked for back to back 
emission and contains the first, second, and higher har- 
monic components. In this respect it is different from the 
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Fig. 4. The correlation function R{A(j>) as function of the relative angle [1]. The solid lines represent the flt using the boosted 
source model with elliptic flow (16), the dashed lines represent correlation for two boosted isotropic sources (13). 



simple — ci cos{A(f)) expression resulting from the momen- 
tum conservation effects. 

If the source is emitting particles in an azimuthally 
asymmetric way in its rest frame 

1 2w2 cos((/)) , (14) 

its distribution in the center of mass frame is 

(1 + 2v2 ((/3± + cos(0))2 - (1 - sin(0)) 

l + /3icos(0) ■ ^ ' 

The correlation function in the relative angle of the pair is 
obtained integrating over the angle of one of the particles 
in the pair 

R{A<j^) = m - 1) l^j g/^(0)/_^(0 + A4>) l) . 

(16) 

The elliptic asymmetry in the particle distribution has the 
same orientation as the momentum exchange (source ve- 
locity). This means, that both the hard scattering, that 
defines the reaction plane, and the generated elliptic flow 
plane have the same orientation. If the generation of the el- 
liptic flow is a collective expansion from fluctuating initial 
conditions [10-13], the orientation of the reaction plane 



from the hard scattering and from the elliptic flow is dif- 
ferent. In that case, the angular distribution in the rest 
frame would be 1 -I- 2v2 cos((/) — A^). An average should 
should be taken over the difference of the two angles A\I/. 
The results are similar, with slightly larger V2, as in the 
scenario using the same orientation of the boost and of 
the elliptic flow (14). 

The formula (16) is fitted to the measured correlation 
function i?(Z\0), adjusting two parameters, the elliptic 
flow coefficient V2, and the boost velocity f3±. The quality 
of the fit is very good (Fig. 4). In Fig. 5 is shown the re- 
sulting elliptic flow coefficient V2- It is smaller than in the 
scenario with momentum conservation effects, studied in 
the previous section. The extracted elliptic flow increases 
with the transverse momentum p±^ and depends weakly on 
the multiplicity. In the p±^ and multiplicity bins where the 
ridge is observed, no significant rise of V2 is seen. The ridge 
appears because the correlations from the transverse boost 
of the source are smaller. In high multiplicity samples the 
boost velocity is much smaller than in the two small mul- 
tiplicity classes (Fig. 6). If such transverse boost velocity 
is present in the particle emitting source, it could mimic 
the transverse flow in the p±^ spectra of particles. 

Our model of background back to back correlations 
from the transverse boost is very simple. We take a sin- 
gle average boost velocity and use massless particles. Go- 
ing beyond these approximations, would require a specific 
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model of the particle emission. The qualitative form of 
the back to back correlation would be similar, but the ex- 
tracted elliptic flow coefficient could depend on the form 
of the dominant angular correlations. 



4 Conclusions 

We show that the angular correlations between particles 
emitted in proton-proton collisions at ^/s = 7 TeV can be 
decomposed into two components. The dominant, back to 
back angular correlations have a kinematic origin. Model- 
ing them with momentum conservation effects results in a 
first harmonic component in the angular distribution. An- 
other scenario assumes a transverse boost of the source, 
this leads to back to back correlations as well, but of a 
different form. The remaining correlations, including the 
small-angle ridge structure, are modeled using an elliptic 
flow contribution. The corresponding elliptic flow coeffi- 
cient V2 is estimated for the first time in proton-proton 
collisions. 

The strength of the elliptic flow depends on the as- 
sumed model for the dominant background, that has to 
be subtracted. In the momentum conservation scenario 
the elliptic flow coefficient is 0.07-0.1 in the bins where 
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Fig. 5. The elliptic flow coeflicient extracted from the CMS 
data using the boosted anisotropic source model. 
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Fig. 6. The transverse velocity of the source extracted from 
two-particle correlation data. 



the ridge is observed, whereas in the transverse boost es- 
timate it is 0.06-0.08. The most conservative estimate is 
obtained using a fit of the second harmonic only in the 
small-angle region of the ridge (lower points in Fig. 2), 
with the result V2 = 0.04-0.06. There are several argu- 
ments showing that the measured elliptic flow contains a 
collective component. It is approximately multiplicity in- 
dependent. The associated structure extends over a long 
range in pseudorapidity. The two particles used to define 
the second cumulant are separated in pseudorapidity, re- 
ducing non-fiow effects. It is tempting to interpret the ob- 
served azimuthally asymmetric fiow as a result of some 
strong rescattering or a hydrodynamic expansion. It re- 
mains as a challenge to provide an alternative microscopic 
explanation, as the available Monte Carlo generators do 
not reproduce the observed structures. For the event gen- 
erators, it would require that such azimuthal correlations 
appear early, between most of the particles in the event. 
We note that there is no additional same-side, ridge-like 
structure in the correlation function. All the observed an- 
gular correlations are perfectly well accounted for using a 
combination of kinematic back to back correlations and of 
the elliptic fiow. 

Finally, we stress again the importance of the possi- 
ble first observation of the elliptic flow in elementary col- 
lisions. It would mean that the short-lived multiparticle 
system created in the collision is very strongly interacting 
and some degree of collectivity appears. 

The author thanks Adam Bzdak and Wei Li for com- 
ments and discussions The work is supported by the Pol- 
ish Ministry of Science and Higher Education grant No. 
N N202 263438. 
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